Mixture of Tonks-Girardeau gas and Fermi gas in one-dimensional optical lattices 
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We study the Bose-Fermi mixture with infinitely boson-boson repulsion and finite boson-Fermion 
repulsion. By using a generalized Jordan- Wigner transformation, we show that the system can be 
mapped to a repulsive Hubbard model and thus can be solved exactly for the case with equal boson 
and fermion masses. By using the Bethe-ansatz solutions, we investigate the ground state properties 
of the mixture system. Our results indicate that the system with commensurate filling n = 1 is a 
charge insulator but still a superfluid with non-vanishing superfluid density. We also briefly discuss 
the case with unequal masses for bosons and fermions. 
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Mixtures of quantum degenerate atoms recently be- 
came a subject of intense studies of both experiment 
and theory. One of particularly interesting systems is 
mixture of ultracold bosonic and fermionic atoms [U-EI, 
which have become accessible through the development 
of sympathetic cooling [l|, Q . The experimental progress 
in manipulating cold atoms in effective one-dimensional 
(ID) waveguides and the ability of tuning the effective 
ID interactions by Feshbach resonance leads experiment 
accessible to the strong correlation regime of ID quantum 
gas 0,01 . Meanwhile, by loading the atomic system into 
the optical lattice @, 0] , one can simulate not only the 
solid state systems in a highly tunable way but also new 
systems which may not be realized in condensed matter, 
such as mixtures of Bose-Fermi atoms. These advances 
open a new channel to investigate numerous phenomena 
of low-dimensional correlated lattice models which play 
important roles in condensed matter physics. 

To gain a deep insight of properties of the low- 
dimensional quantum mixtures, some refined methods 
capable of dealing with strong correlations are especially 
important. For example, the method of Bose-Fermi map- 
ping has been extensively exploited to study the Tonks- 
Girardeau (TG) gas [1, HI . This method has been also 
generalized to study the multi-component quantum gas 
in the infinitely repulsive limit [T3, El- The extended 
Bose-Fermi mapping method is only limited to a spe- 
cial case with no tunable parameter of interaction, where 
all the intra- and inter-component interactions go to in- 
finite. In addition, the Bose-Fermi mixture with equal 
boson-boson and boson-fermion interactions can be ex- 
actly solved by Bethe-ansatz (l2T - [l4l |. Unfortunately, its 
corresponding lattice model is no longer integrable. So 
far, the ground state phase diagram of the ID Bose-Fermi 
Hubbard model in optical lattice has been studie d by 
mean-field theory (THEHi Bosonization method fl7l.ll8j|. 
exact diagonalization method [h| and quantum Monte 
Carlo method [20l - [2^ |. Despite the intensive studies of 



the lattice model [15H201 |22h24| . no anal ytic ally exact 
result has been given except the TG limit [10j, in which 
however the model suffers the problem of a huge degener- 
acy of ground states (GSs) . In this work, we shall study 
the boson-fermion mixtures with the aim to give some ex- 
act conclusions apart from the TG limit and focus on the 
case with an infinite boson-boson repulsion but a tunable 
boson-fermion interaction, which is found to be exactly 
solvable when the hopping amplitude tb for boson equals 
to tf for fermion. 

We consider a mixture system of bosonic and spin- 
polarized fermionic atoms confined in a deep ID optical 
lattice. For sufficiently strong periodic potential and low 
temperatures, the atoms will be confined to the lowest 
Bloch band and the low energy Hamiltonian is described 
by the Hamiltonian 
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where bosonic or fermionic annihilation opera- 

tors localized on site i, and rii a — a. ai a . In principle, 
the interaction parameters U b f and U b can be tuned ex- 
perimentally by the Feshbach resonance. In this work, we 
shall focus on the case with infinitely strong boson-boson 
repulsion, i.e., U b — oo, and a tunable inter-species repul- 
sion Ubf — U. In this case, the boson is a hard-core one 
or a TG gas, for which the states occupied by more than 
one boson are prohibited. Similarly, the states occupied 
by more than one fermion are not permitted due to the 
Pauli principle. However, a boson and a fermion can oc- 
cupy the same site which contributes an on-site energy 
U. In the hard-core limit, the Bose-Fermi mixture model 
can be simplified to 
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Mixture with equal masses. — Firstly, we focus on the 
case that the bosonic and fermionc atoms have the same 
masses, which is approximately satisfied for the Bose- 
Fermi mixture of heavy isotopic atoms, for example, the 
174 Yb- 173 Yb mixture [3[. The Bose- Fermi mixture with 
equal masses provides a solvable limit, which may serve as 
a touchstone for various numerical simulations. For the 
model system with the fermion and the boson having the 
same mass, we have tf, = tf. It is convenient to use the 
following extended Jordan- Wigner ( JW) transformations 
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which maps the Hamiltonian of hard-core Bose and Fermi 
mixture model into a Hubbard model 
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The second mapping in eq. §S§ is introduced to en- 
force the Fermion operators c;^ and fulfil the anti- 
commutation relation {c;^, c^} = 0. The Hamitoni- 
ans of Hbf and Hp have the same spectrum of en- 
ergy. Therefore, we can get the eigen-energy of Hbf 
with tb = tf = t from the well-known Lieb-Wu solution 
of the Hubbard model [25|, i.e., the eigenenergy is given 
by E = —ZtJ^j cos kj with kj determined by the Bethe- 
ansatz equations 
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where 6 n (k) = 2tan- 1 (4fc/n[/), j = !,■■■ ,N, P = 
I,-- - ,M, N = N b + N f , M = N f , N b (N f ) is the 
number of bosons (fermions), and L is the size of the 
system. Here the set {Ij, J a } play the role of quantum 
number. We also solve the ground state energy of Hbf 
by using the numerical exact diagonalization method and 
compare with the result obtained by solving the Bethe- 
ansatz equations. It is found that the numerical result 
agrees with the Bethe-ansatz solution exactly. 

The unitary mapping builds a bridge between the 
hard-core Bose-Fermi mixture and the extensively stud- 
ied Hubbard model. Since these two models sharing the 
same energy spectrum, we can conclude that they have 
the same thermodynamic properties. The GS properties 
of hard-core Bose-Fermi model also share some similari- 
ties with the Hubbard model, i.e., there exists no Mott 
transition from superfluid to Mott insulator for any finite 
U. For the incommensurate filling case, the system is in a 
superfluid phase, whereas the system with the commen- 
surate filling n = 1 is in a Mott phase for any finite U, 
which is characterized by the presence of a charge gap 



and simultaneously a gapless mode of mixture composi- 
tion fluctuations. 

The superfluid density of the bosonic component can 
also be characterized by the bosonic phase stiffness, 
which reflects the response of a superfluid component to 
the imposed phase gradient and is defined as [26| 
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which is proportional to the Drudc weight. Similarly, the 
fermionic stiffness can be represented as 
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Here, Eq is the ground-state energy, (f> b and <pf are the 
component-dependent flux in units oihe/e for the boson 
and fermion respectively, which can be incorporated in 
the wavefunction by making the usual gauge transforma- 
tion <2i i(T — s> e l ^" Til ai }Cr . A finite bosonic or fermionic 
stiffness is characteristic of a superfluid or conductor, 
whereas the stiffness vanishes for a insulator. In the 
presence of component-dependent flux, the eigenenergy 
of the system is also given by E = —It cos kj with kj 
determined by the revised Bethe-ansatz equations 26 1 
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By solving the revised BAEs, we can directly calculate 
the stiffness of system. In order to calculate the charge 
stiffness D c , we set the magnetic flux of hard-core bosons 
(j>b and that of spinless fermions c/)f to be the same, i.e., 
cj)b = 4>f = 4>, whereas the spin stiffness D s is calculated 
by taking (f>b — —4>f — 4>- Explicitly, we have 
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For the case with commensurate filling, we display the 
charge, spin and boson stiffness in (a), (b) and (c) of 
Fig. HI respectively. It is obvious that the charge stiff- 
ness goes to zero quickly when the interaction U exceeds 
a critical value. To extrapolate the critical U c in the 
thermodynamic limit, we make finite size analysis of the 
transition point where the charge stiffness tends to van- 
ish, which is characterized by the minimum of the deriva- 
tive of the charge stiffness as shown in the Fig. Hid). 
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FIG. 1: (color online) The stiffness of system. The panels 
(a), (b), (c) and (d) correspond to the filling factor n = 1, 
whereas the panels (e) and (f) correspond to n = 2/3. Finite 
size analysis of the minimum of the derivative of the charge 
stiffness is shown in inset of (d). Here we have taken Nb = Nf. 



The scaling behaviors of transition points can be fitted 
as U cm = 0.00104+ 1.79789/ InN. When the system-size 
tends to infinity, the critical on-site interactions reads 
U c = 0.00104 ± 0.00354, which covers the zero within 
the scope of fitting error. This means that the system 
is in a Mott phase with zero charge stiffness for any fi- 
nite repulsion between the hard-core boson and fermion 
in the thermodynamic limit. However, as shown in Fig. 
QJb) andQJc), the spin and boson stiffness do not vanish 
even the system is in a Mott phase where the fluctua- 
tion of particle number (charge) on each site is greatly 
suppressed. Furthermore, our results show that the den- 
sity of superfluid fulfills an interesting relation with the 
charge and spin stiffness, Db = (D c + D s )/4 (see ap- 
pendix). The non- vanishing Df, indicates that the boson 
is still in a superfluid phase even in the phase of charge 
insulator. This is induced by the fluctuation of bosons. 
Despite the charge fluctuation is suppressed, the boson 
can tunnel to the neighboring sites by a virtual second 
order process. For the incommensurate case with n < 1, 
both the charge and spin stiffness do not vanish as shown 
in Fig. [lje) andQJf), and no a Mott phase exists. 

Particular attention should be paid to the special cases 
with a partial commensurate filling, for example, cases 
with rif = 1, n;, < 1. Then the system favors that every 
site has one fermion which supplies a background and the 



bosons can hop on it freely. The energy of this state is 
tL f kb 

E = N b U - — / cos(fc + 4> b /L)dk, 
77 J-k b 

where fcf, = 7rn&. It is straightforward to get Db = 
t sm(Tmb) / tt and Df — 0, which indicates that the bosons 
form a superfluid whereas the fermions form an insulator. 
Similarly, for the case with rib = 1 and n./ < 1, the bosons 
are in an insulator state which provides a homogeneous 
background for the fermions which form a conductor with 
Df = t sm(Tm f) / n . 

Despite Hbf and Hp sharing the same energy level 
structure, they have different ground state wavefunction 
due to the intrinsically different exchange symmetry of 
the wave functions for Bose and Fermi systems. Sup- 
posed the wavefunction of the Fermi Hubbard model 
is given by ^f(xi, •• ■ , x n ; x n +\, ■ ■ ■ , x at), the wavefunc- 
tion of Hhb can be constructed as 

^bf(xi,--- , x n ;y n +i, ■ ■ ■ ,vn) 
= Y\_ s 9 n ( Xi _ x j)^F{xi, ■ ■ ■ ,x n ;y n+1 , ■ ■ ■ ,y N ), 

i<j 

where sgn(Xi — Xj) — (x, — Xj)/\xi — Xj\ is the sign 
function. Consequently, the observable associated with 
the wave functions rather than the energy level struc- 
tures should display different behaviors, which can be 
displayed in the off-diagonal density matrix and the mo- 
mentum distributions of the hard-core boson. Explicitly, 
the density matrices of boson and fermion are defined 

as pfj = (al >b a,jfi\ and pfj = (aljdjj) respectively, 

which exhibit quite different behaviors for boson and 
fermion. The momentum distribution can be obtained 
by the Fourier transformation of the corresponding den- 
sity matrix. For example, we have 

In Fig. ([%]). we display the momentum distribution of the 
Bose and Fermi systems respectively for the case with 
commensurate filling rib — n f — 1/2- It is shown that 
the momentum distribution of the hard-core bosons has 
a sharp peak, which reflects the bosonic nature of the 
particles. With the increase in inter-component inter- 
actions, the momentum distribution spreads wider and 
wider but the pronounced peak around the zero momen- 
tum is kept. While the momentum distribution of the 
fermion exhibits the free Fermi distribution for U = 0, 
it becomes more wider and develops a wide tail with the 
increase in U . 

Mixture with unequal masses. — Finally, we give a brief 
discussion on the case with % ^ tf which corresponds to 
the system where the bosonic and fermionc atoms have 
different masses, such as mixture of 7 Li and 40 K. In gen- 
eral, the single particle hopping amplitude is inversely 
proportional to the mass of atoms, i.e., tf/tb = nrib/mf. 
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FIG. 2: (color online) Momentum distributions n B (k) for 
hard-core bosons (a) and n F (k) for spinless fermions (b). Here 
L = 14, iV b = 7 and 7V f = 7. 
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FIG. 3: Phase diagram for the asymmetric Bose-Fermi mix- 
ture with n = 4/5. 



Taking the mixture 7 Li and 40 K as an example, we have 
tp = tt/tb ~ 0.175. For the general case with tp ^ 1, the 
mixture of hard-core bosons and fermions is no longer 
exactly solvable, although a mapping to the asymmetric 
Hubbard model [27j still holds true via the generalized 
JW transformation. 

In the heavy Fermi mass limit tp — >■ 0, the fermions lose 
the mobility and the asymmetric mixture model is related 
to the Falicov-Kimball model [28] via the generalized JW 
transformation. Away from integer-filling with n < 1 , the 
heavy fermion and the light bosons are favorable to stay 
in different regimes to lower the kinetic energy of bosons. 
Therefore the system is expected to have a transition 
from density wave to phase separation. Next we shall 
use the density matrix renormalization group method to 
obtain a quantitative phase diagram for the mixture of 
TG gas and Fermi gas with Nb = Nt and a filling n = 
4/5. 

Taking into account that the dominating configura- 
tion of fermionic atoms is quite different in the density- 
wave phase or phase with phase separation, we introduce 
the following structure factor of density wave (DW) of 



fermionic atoms 



(10) 



where q = 2nir/L and n = 0, 1, ••• ,L. We calcu- 
late the structure factor as a function of tp for different 
modes for systems with different sizes L and strengthes 
of interactions U. The results show an obvious com- 
petition between the modes of Sfdw(<Z = 2ir/L) and 
Sfdw(q = Nit/L). In the heavy fermionic atom limit 
with small tp, <Sfdw(<7 = 27r/L) dominates, which indi- 
cates phase separation in this region [29| where configu- 
rations of fermionic atoms like ]/, /, /, /, o, o, o, o, o, o) are 
found to be dominant. On the other hand, as tp 1, 
Sfdw(<1 — Nir/L) exceeds Sfdw(<Z = 27r/L), which im- 
plies that fermionic atoms distribute uniformly on the 
optical lattice. Then together with bosonic atoms, the 
ground state becomes the so called DW state similar to 
the state of symmetric model, which is the limiting case 
with tp = 1 studied in the above section. Consequently, 
we can determine the transition point on the U — tp plane 
for a finite system from the intersection of the structure 
factor of two modes. In Fig. [31 we plot the phase dia- 
gram on the U — tp plane for systems with a filling factor 
n = 4/5 for different system-sizes L = 15,20,25. The 
infinite size limit is obtained by extrapolation from the 
finite-size analysis. Below the phase boundary line, the 
phase-separation phase dominates, while the DW phase 
is dominated above the boundary line. 

In summary, we study the mixture of TG gas and 
fermions in a ID optical lattice. We show that this sys- 
tem can be mapped to the Fermi Hubbard model by 
a generalized Jordan- Wigner transformation and thus is 
exactly solvable when the bosons and fermions have the 
same masses. Based on the Bethe-ansatz solution, we 
calculate the charge stiffness and the superfluid density 
for the systems with either commensurate or incommen- 
surate filling. Our results show that the system with 
commensurate filling factor n — 1 is a charge insulator 
but remains to be a superfluid characterized by a non- 
vanishing superfluid density. We also give a brief discus- 
sion to the case with unequal boson and fermion masses. 
In the heavy Fermi mass limit, our result indicates that a 
phase-separation phase arise. The phase transition from 
the density wave phase to the phase-separation phase is 
discussed. 
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Appendix A 

In this appendix, we provide a derivation of the re- 
lation D\, = \(D C + D s ). In the presence of flux, the 
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Hamiltonian (2) becomes 

H BF = -t^ (4,a i+ i 6 e^>/ L +H.c) 

i 

i i 

Expanding the Hamiltonian in terms of <p b /L and 4>f /L, 
we have 

+Uj2n i , b m,f + 0(cPi,cj ) j), 

i 

where T a = -tJ2i {a\a a i+i° + H - c -) and 3a = 

{ a la a i+^ ~ H.c.^j with a = b,f. So according 
the perturbation theory, 
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where the current operators: j c = j b + jf, j s = j b - jf, 
and T C = T S = T } +T b = 2T b . Then 
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Finally, we have 
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D b = -(D C + D s ). 
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